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Abstract 

We study the three-dimensional generalized six-state clock model at val- 
ues of the energy parameters, at which the system is considered to have the 
same behavior as the stacked triangular antiferromagnetic Ising model and 
the three-state antiferromagnetic Potts model. First, we investigate ordered 
phases by using the Monte Carlo twist method (MCTM). We confirmed the 
existence of an incompletely ordered phase (lOPl) at intermediate temper- 
ature, besides the completely ordered phase (COP) at low-temperature. In 
this intermediate phase, two neighboring states of the six-state model mix, 
while one of them is selected in the low temperature phase. We examine the 
fluctuation the mixing rate of the two states in lOPl and clarify that the 
mixing rate is very stable around 1:1. 

The high temperature phase transition is investigated by using non- 
equilibrium relaxation method (NERM). We estimate the critical exponents 
[3 = 0.34(1) and u = 0.66(4). These values are consistent with the 3D-XY 
universality class. The low temperature phase transition is found to be of 
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first-order by using MCTM and the finite-size-scaling analysis. 

I. INTRODUCTION 

The systems of Zg symmetry in the three dimensions (3D) have attracted interests be- 
cause of their peculiar nature of orderings. In particular, there have been many works on 
the nature of intermediate temperature phases. 

In two dimensions, the six state clock model shows three phases. At intermediate tem- 
perature, the discreteness of the six states is irrelevant and the system shows the same 
properties as that of XY model, which is the Kosterlitz-Thouless phase where the long 
range order does not exist but the correlation length of the order diverges 

In three dimensions, the six state clock model also shown natures which are similar to the 
ordered phase of the XY model. However it found that there exists only a phase where one of 
the six state of the model is chosen as the long range order 0. The apparent behavior is due 
to the large fluctuation at the higher temperature region of the ordered phase. However, if 
we generalize the energy structure of the model, a new type of intermediate phases appears. 

Ueno et al have introduced three-dimensional generalized six-state clock model (3D- 
6GCL model) 0. The 3D-6GCL model can be regarded as prototypes of wide range of 
Zg models. That is, it represents various categories of the Zg models according to the 
values of the energy parameters of the model. The existence of intermediate phase have 
attracted interest for the following models: three-dimensional three-state antiferromagnetic 
Potts model (3D-3AFP model) |Q, and stacked triangular antiferromagnetic Ising model 
(STAFI model) 0. Although, the existence of various type of intermediate phases for 
these models have been proposed, there is no general understanding on the nature of the 
intermediate phase(s) p|-p!2|. If we consider only one universality class for models with Zg 
symmetry, we conclude that the intermediate phase is an apparent phase as we found in the 
regular six state model. However, it has been clarified that the situation is not so simple and 
interesting new type of phase exists in the generalized case. In fact, the 3D-6GCL model was 
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introduced in order to understand general properties of other Zg models [Q. In the present 
paper, we study the 3D-6GCL model with a set of the parameters which is considered to 
correspond to the STAFI model and the 3D-3AFP model. 
The Hamiltonian of the 3D-6GCL model is given by 

(1) 

ihj) 

where runs over nearest neighbor pairs, and Ui is the spin variable which takes one 
of 1,2, ■■■ and 6. We set Sq = Q > ei > €2 > as in Fig.0. Namely, we only study 
ferromagnet cases where all the spins occupy the same state in the ground state. We expect 
that this ferromagnetic ordered phase appears at the low temperature. This ordered phase is 
called completely ordered phase (COP). Since the COP has six degenerate states, the order 
parameter space of this model is illustrated by a hexagon as in Fig.Q. The bold solid line 
represents the one of the six COP states. The dotted line represents a mixture of neighboring 
two COP states which is a candidate of an intermediate phase, which we call incompletely 
ordered phase (lOPl). As to the distribution of the spin state in the incompletely ordered 
phase, there have been proposed the two kinds depicted in Fig.|^. The lOPl discussed 
above corresponds to Fig.^(a). The type of Fig.^(b) was called I0P2 which was one of the 
possible mixed state. For Si = 0, the COP does not appear because the ground state is 
macroscopically degenerate. This degenerate state is lOPl and corresponds to the ground 
state of the 3AFP model. Taking the approach of physical percolation to phase transitions, 
Ueno rigorously proved that for si = 0, at least one kind of lOPs can exist at T < £2 = ^3 
pr3| . For El > , Ueno and Kasono proposed that this model has two incompletely ordered 



phases (lOPl, I0P2) besides the COP by using the Monte Carlo twist method (MCTM) 
In this paper, we reexamine the nature of the intermediate phase(s) in detail. 
On the 3D-3AFP model, so far various types of phases at intermediate temperatures have 
been proposed in the previous studies, such as the permutationally symmetric sublattice 
phase (PSS phase) which corresponds the I0P2 or the gapless phase which is similar to 
the ordered phase of the three dimensional XY model which we call XY phase [^. Recent 



theoretical and numerical studies on the 3D-3AFP model have revealed that the existence of 



these phases are fake owing to the finite size effect |H|-[T6[] . There is only a low-temperature 
phase that corresponds to the lOPl. 

We will show that the correlation length of the fluctuation of the order parameter in 
lOPl is finite. However the correlation length just below the critical point is so long that 
the phase looks like the ordered phase of the 3D XY model where the correlation length is 
infinite although the long range order exists. In order to distinguish the difference between 
the XY phase and lOPl we have to be careful in the examination of properties of the phase. 

In this paper, we study a 3D-6GCL model by using the Monte Carlo (MC) simulation 
and clarify properties of the phase at intermediate temperature and natures of the phase 
transitions. In this study, we choose the energy parameters of the 3D-6GCL model as eg = 0, 
£1 = 0.1 and S2 = s^ = 1.0 which are the same as those in IQ. 

In order to examine the possibility of many ordered phases by the MCTM, one needs 
a full set of appropriate boundary conditions which cause different types of domain-walls, 
and calculate size dependence of the excess free energy. In the present study, we deny the 
existence of 10 P2 as intermediate phase from the dependence of the excess free energy on 
the boundary conditions. The intermediate phase shows apparent property of the XY phase 
where Zg discreteness is irrelevant. However, we also deny the existence of XY phase in the 
thermodynamic limit. This apparent property is attributed to the long correlation length 
at the high-temperature region of the intermediate phase. Finally, we conclude that only 
lOPl exists as the intermediate temperature. Appearance of the fake XY phase near the 
critical point has been also found in the six-state clock model, where the Zg discreteness is 
very weak in the ordered phase near the critical point. However the Zg discreteness is always 
relevant in the ordered phase 0. 

Moreover, we examine the properties of the phase transitions. We consider that the 
high temperature phase transition belongs to 3D-XY universality class. In order to confirm 
whether this transition belongs to 3D-XY universality class, we study this transition by 
using non-equilibrium relaxation method (NERM) and estimate the critical temperature 



and critical exponents. Ueno and Kasono pointed out that the low-temperature transition 
between the lOPl and the COP is of the first order because the two phases have no symmetry 
relation of group and subgroup [^]. We investigate this phase transition by using MCTM 
and the usual finite size scaling analysis, and confirmed that it is of the first order. 

The outline of this paper is as follows. In Sec. II, we study properties of the intermediate 
phase. In Sec. Ill, we study that the high temperature transition and in Sec. IV, we investigate 
the low temperature phase transition, and the Sec.V is for summary and discussion. 



II. INTERMEDIATE PHASE 

The properties of the phases of the 3D-6GCL model have been studied by the MCTM [Q, 
which is a simulation method under special boundary conditions to detect the domain-wall 
excess free energy. We prepare systems with fixed boundaries in one direction and impose a 
periodic boundary condition in the other directions. As a reference system, we set the same 
state a at the boundaries (a-a) in which we expect that the ordered phase a appears. We 
also prepare a system where we set the states a and (3 at the boundaries in which we 

expect an interface. The interfacial free energy is defined as the excess free energy between 
the two systems: 

AF"^ (T, L) = F"^ (T, L) - F"" (T, L) , (2) 

where L is a liner size of the system. Here we consider an L x L x L system. The size 
dependence of AF°'^(T, L) is given by the following asymptotic behavior for large L 

AF(T, L) ~ A{a, p)L^'^^^^\ (3) 

where ipapiT) is called stiffness exponent. The ordered phase is classified by the value of 
stiffness exponent as 

D — 1 for a domain- wall type interface < oo) 
i^afsiT) = \ D — 2 for a gepless interface (spin- wave type ^ = oo) • (4) 
non-integer for a new type on interface 



When the interface does not appeare, the stiffness exponent becomes negative because 
AF(T, L) for L ^ oo. 

In the present work, we use three sets of boundary conditions which are shown in Fig.^. 
There are two kinds of boundary conditions; in 0i we fix the boundaries to be two of the 
low temperature phases. The relative angle (twist angle) between the two phases can be 
^, Y' ^- Here we take the cases | and n which are shown in Fig.^ (a) and ^ (b), 
respectively. In 02 we fix the boundaries to be two of the lOPls in Fig.^ (c) (twist angle 
7r/3). As shown in Table H, we distinguish the five phases from the signs and the values of 
the stiffness exponent. We perform MC simulation for the lattices with liner size L = 8, 10, 
and 12 with boundary conditions 0i(vr), 0i(7r/3) and 02(7r/3). At each temperature, after 
discarding first 5000L, 5000L and lOOOOL steps, we calculate quantities of interest using 
data of next 25000L, 25000L and 50000L steps for the three sets of boundary conditions, 
respectively. 

In Fig.^, we show the temperature dependence of ip under the boundary conditions. From 
this figure, we may conclude the existence of the disorder phase, the XY phase, the lOPl and 
the COP, as the temperature decreases. The high temperature transition point is estimated 
to be Tq ~ 1.35(5), and the low temperature transition point is to be Tp ~ 0.30(5). We see 
a change around T ~ 0.7 which can be a phase transition between the XY phase and lOPl. 

In order to confirm the region of phases, we investigate the size dependence of the 
stiffness exponents for 0i(7r/3) for larger sizes. The region of the XY phase decreases as 
L becomes large (Fig.^. The size- dependence strongly suggests that the region of the 
XY phase disappears in the limit L — oo and only the lOPl exists in the intermediate 
temperature region. 

If we consider the domain wall between the ordered states of the type of lOPl, there is 
no domain wall in the (pilir/S). On the other hand, one domain wall appears in the case of 
02(7r/3), and two domain walls appear in the 0i(7r). Therefore, stiffness exponent for the 
cases 01 (tt) and 02(7r/3) should be the same. However, we notice that the difference of the 
values of stiffness exponents of 0i(7r) and 02(vr/3) is significant at the lOPl as we see at 

6 



Fig.^. The reason for this difference is that the effect of twisting in the system by 0i(7r) is 
stronger than one of the twisting by (/)2(7r/3). Actually, when we examine ipiiii) and ■?/'2(7r/3) 
in much larger lattice sizes and we confirm that the values of the stiffness exponents ^/'i(vr) 
and ^/'2(7i"/3) approach the same value ip = 2 which indicates the domain- wall type order. 
In Fig.^ we show a snapshot of the domain- wall of the boundary in 02(7i"/3) of L = 160 at 
T = 0.4, where we find a localized domain-wall clearly. 

Next, following the idea of Oshikawa |]14| for the 3D-3AFP model, we make the finite 
size scaling plots of — (cos6^^) by ordinary Monte Carlo method on lattices L = 12, 18,24 
and 30. Here, —(cos 66') is the order parameter for the Zg symmetry breaking and 6 is the 
angle of magnetization from one of the COP states. After discarding first lOOOOL steps, we 
calculate quantities of interest for next 40000L steps at each temperature. In Fig.|^, we show 
the scaling plots in the form of the scaling function given by Oshikawa 

{cosQe)^f{cL'{Tc-Tr\y'^\). (5) 

Here, we put the value of the critical temperature Tq = 1.384 which was obtained by using 
the NERM which will be explain a later. We obtain a good scaling plot when we use 
^luel = 3.2. It is important to notice here that behavior of the order in the region of the 
lOPl is expressed by a scaling function using the high temperature critical point Tq. This 
fact indicates that other phases do not exist between the lOPl and the disordered phase 
because the fluctuation in the lOPl phase is controlled by the high temperature critical 
point. However, it is very difficult to estimate the accurate value of this exponent because 
the strong finite size effect appears in the intermediate phase even at the largest system 
which we can calculate. After all, we conclude that the 3D-6GCL model with Eq = 0, 
El = 0.1 and E2 = Es = 1.0 has two phases, i.e. lOPl as the intermediate phase and COP 
as the low temperature phase. 

In order to investigate the nature of lOPl, we examine the following quantity; 

X = ^MiPi - P2r) - {{pi - P2))% (6) 
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where pi denotes the number of the spin variables which take the value of the i-th state. 
We adopt the boundary conditions 02(0). The asymptotical form of x is written x ~ AL^. 
The exponent x depends on the type of the order. When the fluctuation is of the XY phase, 
we expect x = 2. On the other hand, when the fluctuation has a finite correlation length, 
x = 0. 

We perform MC simulations for the lattices with the liner sizes L = 8, 10, ■ ■ ■ and 24. 
For each temperature, after discarding first lOOOOL steps, we calculate x using data of next 
50000L steps. In Fig.^, we show the size and temperature dependence of x. We find that 
the exponent x decreases in the lOPl region. Moreover, we perform MC simulation for the 



lattice up to the liner size L = 80 at T = 0.4 (in the lOPl region). In Fig.|TO|, we show 
the temperature dependence of x at T = 0.4. The value of x approaches the constant as L 
increases. Thus we consider that the value of x decreases and approaches as L — oo, and 
we conclude that the fluctuation in lOPl is the type of the domain- wall and the correlation 
length of the fluctuation of the order parameter is finite. 

III. HIGH TEMPERATURE PHASE TRANSITION 

The high temperature phase transition in Zg models belongs to the 3D-XY universality 



class Indeed, the transitions in the 3D-3AFP model and the six-state clock model 
are found to belong to the 3D-XY universality class. We consider that the high temperature 
phase transition in the 3D-6GCL model also belongs to 3D-XY universality class. 



We investigate this phase transition by the NERM p7|-[T9|. The NERM is an efficient 
numerical technique to estimate the critical point and critical exponents from a dynamical 
process toward the equilibrium state from an ordered state. The decay of the order pa- 
rameter m{t) shows a power-low only at the critical point. Detecting such point, a precise 
determination of the critical temperature can be done. The asymptotical form of m{t) at 
the critical temperature is written as 

m{t)r^t-^^. (7) 



To examine asymptotic behavior of m{t) clearly, we introduce a local exponent Xm{t) 

alogt 

This exponent Xm{t) corresponds to f3/zv. Further, we consider the following function and 
the local exponent. 



^ d\0gfmm{t) . . 

alogt 



where \mrn{t) corresponds to d/z. Therefore, we obtain the exponents 13/v and z indepen- 
dently from these quantities. We simulate the relaxation process starting from an initial 
state which is set to be the lOPl state, and measure the magnetization m{t). In Fig.|Tl|, 
we show Am(t) for the temperatures near the Tq. Here, MC simulation is performed in a 
lattice L = 60. We show the MC steps and the number of samples in Table The curve 
for T = 1.385 turns up while the curve for T = 1.383 turns down. Therefore we locate the 
transition temperature in 1.383 < Tq < 1.385, denoting Tq = 1.384(1). 

At this accurate value of Tq, we calculate the relaxation of quantities. We perform about 
160000 independent runs up to 200MCS and average the process to obtain m(t). We obtain 
the critical exponents P/u and z by making use of (|) and (|^). In Figs.[T^ and |IB|, we show 
z(t) and P/uit), respectively. From the extrapolated values of z(t) and P/uit), to 1/t = 0, 
we estimate z = 2.05(5) and p/u = 0.515(10). 

Moreover, assuming Tq = 1.384(1) and P/v = 0.515(10), we estimate P from the scaling 
plots of data of the magnetization which are obtained from the ordinary equilibrium MC 
simulation. Seeking the value of P by which the date collapse into a scaling function, we 
estimate p = 0.34(1) for the best fit (Fig.|14D. Since the values, P = 0.34(1) and u = 0.66(4), 
are close to the 3D-XY universality class {P = 0.345, u = 0.669), we conclude that this 
phase transition belongs to the 3D-XY universality class. 
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IV. LOW TEMPERATURE PHASE TRANSITION 



While many researchers have studied the high temperature transition in the Zg model, 
but few study has been done for the low temperature one. Ueno and Kasono pointed out 
that the low temperature transition between the lOPl and the COP is of first-order because 
they have no symmetry relation between group and subgroup which is usually seen in the 
second-order transition. 

In this paper we study this problem by a finite-size-scaling analysis of (cos6^^). If the 
transition is of first order, the order parameter is scaled by (cos 60) ~ /(eL*^) where e = 
(T — Tq)/Tq, d is the dimension of the system and f{x) is a finite size scaling function. We 



obtain a good scaling plot at Tp = 0.295(5) depicted Fig.[T^. Thus we confirm this phase 
transition is of first-order. 



We also confirm this result by an analysis of extension of the MCTM p0|j21[] . Now, we 
study the interfacial energy AE is defined in the same way for the excess energy due to 
the boundary condition causing the interface, eq.(2). In a completely ordered phase, the 
size dependence of the excess energy becomes AE ~ L"^^^. At the critical point, it becomes 
AE ~ AL^/'^ . We define the stiffness exponent of the energy ipE as AE ~ AL'^'^ . For the 3D 
Ising model at the critical point ipE = 1-59 which is considered to be a large value among the 
3D models exhibiting the second order phase transition. We study ipE for the present phase 
transition. We use the boundary condition 0i(7r/3). In Fig.|l6|, we show the temperature 
and the size dependence of ip-^. We confirm that ip-^ approach to 3 which is much larger 
than the values oi ipE oi the second order phase transition in 3D models. For the first order 



phase transition, it is known that ipE = D from the argument of wetting [^. This result 
also suggests that the present phase transition is of first order. 
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V. SUMMARY AND DISCUSSION 



We study successive phase transitions in the 3D-6GCL model and found an intermediate 
phase and two phase transitions. First, it has been revealed that the intermediate phase 
is single phase of the lOPl. In the high temperature region of the intermediate phase, the 
correlation length is very large and the system shows an apparent XY behavior. In this 
region the stiffness constant is very sensitive to the boundary condition. Therefore, the 
wrong results were obtained in the previous numerical studies Q . The present results agree 
with the phase diagram of the 3D-3AFP model obtained by Kishi and Ueno, and our results 



of the STAFI model |2^. Thus, we have obtained the same intermediate phase in the three 
3D Zg models, i.e. the 3D-6GCL model, the 3D-3AFP model and STAFI model. Moreover, 
we found that the mixing rate of the two states in the intermediate phase is steady and 
does not show anomaly fluctuation of the rate proposed preciously 0. This lOPl with the 
steady mixing of two COP states can be regarded as partially disordered phase proposed by 
Mekata in the original work on the triangular antiferromagnets. Here we understand 
that the entropy effect due to the flustration allows different type of intermediate phases in 
two and three dimensional model. 

Secondly, we examine properties of the high temperature phase transition, and the tran- 
sition temperature is estimated to be Tq = 1.384(1) with the critical exponents f3 = 0.34(1), 
u = 0.66(4) and z = 2.05(5). These values of the exponents are close to those of 3D-XY 
universality class. We conclude that the universality class of this transition belongs to the 
3D-XY universality class. 

Thirdly, we investigate the low temperature phase transition in the 3D-6GCL model by 
using the finite-size-scaling analysis and the MCTM. All of the results support that this 
transition is of first-order. This first-order phase transition is a transition at which the 
(cos 60) changes from —1 to 1 discontinuity. 
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TABLES 



TABLE I. Values and signs of stiffness exponents for each phases 





Disorder 


lOPl 


I0P2 


XY 


COP 


(/)i(7r) 




+ 


+ 


1 


2 


<^i(7r/3) 






+ 


1 


2 


02(7r/3) 




+ 




1 





TABLE II. Details of the simulations of NERM for the 3D-6GCL model. 



lattice size 


temperature 


Monte Carlo step 


the number of samples 


60 


1.383 


200 


20000 




1.384 


200 


60000 




1.385 


200 


60000 




1.386 


200 


20000 
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FIGURES 



FIG. 1. Energy level of a neighboring spin pair of the 6GCL model 

FIG. 2. Order parameter space in 3D-6GCL model. The solid line connecting a vertex and the 
center is one of six low temperature state and the dotted line between the neighboring solid lines 
is one of six lOPl states. 

FIG. 3. Schematic of the one-spin distribution functions for two kinds of lOP's (a) lOPl and 
(b) I0P2. 

FIG. 4. Three boundary conditions for investigation of the intermediate phase by using the 
MCTM. (a) 0i(7r), (b) 0i(7r/3) and (c) (/)2(7r/3). 

FIG. 5. Temperature dependence of the under the three boundary conditions. 

FIG. 6. Temperature and size dependence of the stiffness exponent under the boundary condi- 
tion ^i(7r/3). 

FIG. 7. Scaling plots of —(cos 6^) for various system sizes and with iy\yQ\ = 3.2 and c = 0.065, 
and we use Tq = 1.384 that is estimated by using the NERM. is the modified Bessel function. 

FIG. 8. Snapshot of one of layers under the boundary condition (/)2(7r/3) at T = 0.4. The colors 
represent the states. The green cells is the state 1, The red cells is the state 2, and The blue cells 
is the state 3. There are some of cells of different states by fluctuation. 

FIG. 9. Temperature and size dependence of x 

FIG. 10. Size dependence of x at T = 0.4 

FIG. 11. Relaxation of Xm{t) 
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FIG. 12. Relaxation of the z{t) at Tq = 1.384 

FIG. 13. Relaxation of the P/i^{t) at Tq = 1.384 

FIG. 14. Scaling plots of m for various system size with /3 = 0.34, and we use Tq = 1.384 and 
(3 /u = 0.515 that is estimated by using the NERM. 

FIG. 15. Scaling plots of (cos 66*) around the low temperature phase transition with Tp = 0.295 
FIG. 16. Temperature and size dependence of the stiffness exponent of AE 



17 



